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Abstract 

We study the three-body baryonic B decays of B ^ pp{K*,p) in the standard model. The 
baryonic matrix elements are calculated in terms of the SU{3) flavor symmetry and the QCD 
power counting rules within the the perturbative QCD. We flnd that the decay branching ratios, 
angular and direct CP asymmetries of {B~ — > ppK*~ , B^ ppK*^ , B~ — > ppp~) are around 
(6, 1, 30) X 10~®, (13, —27, 11)% and (22, 1, —3)%, which are consistent with the current BaBar and 
Belle data, respectively. The large values of the branching ratio in B~ ^ PPP~ and the direct CP 
asymmetry in B^ —>■ ppK*^ are useful to test the standard model and search for new physics. 

PACS numbers: 13.25.Hw, 12.38.Bx, ll.30.Er, 13.87.Fh 
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I. INTRODUCTION 



The three-body baryonic B decays of ApvT"*", B~ Ap7, B ppK, B~ ppTT~ 

and B^ h.KK~ have been observed by BaBar and Belle Collaborations at the levels of 
10"^ on the decay branching ratios with near threshold enhancements in the 

dibaryon invariant mass spectra. On the other hand, it is still inconclusive whether similar 
threshold enhancements exist in the modes with the vector mesons such as i? — > ppK* due 
to the large experimental uncertainties in the present data, given by 

m 



Br{B ^ ppK 



[lO.Stlitl-',) X 10- 



Br{B^ 



(Belle) , 

= (4.94 ± 1.66 ± 1.00) X 10^^ (Babar) 
ppK*^) < 7.6 X 10-^ 90% C.L. (Belle) , 

= (1.28±0.56tJ]:}?) X 10"^ (Babar). 



(1) 



Theoretically, the dibaryon threshold enhancements in the three-body baryonic B decays 
were first conjectured in Ref. [8|. Subsequently, various interpretations, including models 
with a baryon-antibaryon bound state or baryonium 9|, exotic glueball states 10, frag- 



mentation 



ll| and final state interactions 12] have been proposed. Nevertheless, only the 
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15] and the QCD counting rules within the framework of 



18| lead to consistent and 
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systematic 
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21 



understandings 



approachs of the pole model [13| 

the perturbative QCD (PQCD) ^ ^ 

on B ^ BB'M (B, B' = A and M = D,K, vr, 7) Q, Q; 
However, for the vector meson mode of B^ ppK*~ , the pole model predicts that a branch- 
ing ratio smaller than Br{B~ — ^ ppK^), which does not agree with the data as pointed out 



22 



23 



24 



25]. 



in Ref. [15|, while the corresponding PQCD study has not been done in the literature yet. 
In this paper, we give a systematic analysis on B ^ ppMy with My = K* and p based on 
the B pp transition hadronic matrix elements parametrized in Ref. 24]. 

The paper is organized as followed. In Sec. II, we present the formulation for the 
decays oi B ppMy in the standard model. We first derive the decay amplitudes via the 
effective Hamiltonian and parametrize the hadronic matrix elements in terms of various form 
factors. We then list the formulas for the decay widths as well as the angular and direct CP 
asymmetries. In Sec. Ill, we show the numerical analysis and summarize our results in Sec. 
IV. 
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II. FORMULATION 



From the effective Hamiltonian at the quaA level 26|, |27|, the decay amplitude of B 
ppMy is separated into two parts, given by 

A{B ppMv) = C{B ppMv) + T{B ppMy) , 



C{B ^ ppM^ 



T{B ppMv) 



G p 

7! 



yubV*qa2{pp\{uu)y_A\^) 



ytbVtq[ a^{pp\{uu 



dd + ss)v-a\^) 



-ai{pp\{qq)v~A\^) + a^{pp\{uu + dd + ss)v+a\Q) 



~{pp\{2uu — dd — ss)v-a\^) 



{My\{qh)v^A\B) 



+VtbV:^2a, {pp\ {qq) 5+P 1 0) {My \ms~p\B) 

B), 



G 

-^aMy {Mv\{qq')v-A\0){pp\ (q'b) v-A 



72 



(2) 



where Gp is the Fermi constant, Vq.q. are the CKM matrix elements, {qiqj)v±A = ± 
75)g„ {m)s±P = q^{l ± 75)g, and a. = cf ^ + ctH/N, with the color number iV, for z =odd 



26 



23]. In Eq. 



(even) in terms of the effective Wilson coefficients cl^'^^ , defined in Refs. 
q = s{d) and q' = u{d) correspond to i? ^ ppK* (p) and the charged (neutral) modes, and 
c^Mv with My = K*~ , p~ , K*^ and p° are given by 



aK*- 
otp- 



VubV*aa2 + VtbVta { a4 - -ag 



(3) 



respectively. We note that we have used the generalized factorization method with the 
non-factorizable effects absorbed in N^^^. 

To evaluate the amplitude in Eq. ([2]), we need to know the B My transition form 
factors and the vector-meson decay constants as well as the time-like baryonic and B — >■ 
pp transition form factors due to the vector, axial-vector, scalar and pseudoscalar quark 
currents, denoted as V^, A^, S and P, respectively. 
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The mesonic B — > My transitions have the following structures: 



2Vi 



rriB + 



{Mv\A^\B) 



6 ■ q 

{niB + mMy)Ai + i—^q^{2mMv)AQ 



{Pb + Pa/ 



— m 



ie* ■ q) 



A, 



niB + 



(Mv-lVjO) = fMynxMySl 



(4) 



(5) 



where q = Ps — PMv^ ^-iid Ao,i,2 are the form factors of the B — > My transition, and e* 
and /mv- are the polarization and the decay constant of the vector meson, respectively. 
For the time-like baryonic form factors, we have 

F2 



{pp\v,\o) = u!^Fa,+ 



nip + nip 



icr^^uiPp + Pp)fi }v 

F2 



nip + nip 



{Pp-Pp)t, >v , 



{pp\A^,\0) = u^gAif. 
{pp\S\0) = fsuv , 
{pp\P\0) = gpu-f^v, 



hA 



nip + nij 



(6) 



where u and v are the spinors of the baryon pair, while Fi, F2, qa, hA, fs and qp are the 



form factors. Based on the QCD counting rules in the PQCD approach [16, 
can be parameterized as 

1 „ „ 



17 



-^1 = 2<^|| + 2% , 9a = -G\\ 



:Gtt , gp — fs 



for {pp\{uu)v,A,s,p\^) and 



2l|,they 



(7) 



Fi = -Gil + -Gtv, qa = -Gil 

3 II 3 II ' 3 II 



g^ij ) gp — fs 



for {pp\{dd)Y^A,s,p\^) ^ where G\\, G-jj and -F|| are the functions of t expanded by 



G 



t2 



t 



1 -7 



Gn 



5 

t2 



t 



t2 



t 



1 ~7 



(8) 



(9) 



with 7 = 2.148, Aq = 300 MeV, where the power expansion of 1/t^ is due to 2 gluon 
propagators attaching to valence quarks to form a baryon pair. In Eq. ([6]), F2 is suppressed 
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by l/(tln[t/AQ]) in comparison with Fi 28|, |29|] and therefore can be safely ignored; while 
Ha is given by 



IT-A = gA 



(10) 



by equation of motion. 



For B ^ pp transition form factors, the general forms are given by 241 



{pp\V^\B) = iulgij^ + g2i(Tf,uP'' + gsPf, + g^Pp + Pp)t^ + g^{Pp - Pp)a.]75W , 
{pp\A^\B) = m[/i7/, + f2i(Tf,uP'' + /sPm + MPp + Pp)i^ + MPp - , (H) 

where p = Pb — Pp — P p is the emitted four momentum. By using the SU(3) flavor and 
SU(2) spin symmetries 20|, |2J], they are composed of another set of form factors in the 



QCD counting rules with the relations of 



gi 



1 , 

2%) Ji 



for {pp\{ub)v,A\B ) and 



gfi'll ) gj — ^f\\ — fj 



-n = -h 



(12) 



gi = 3^11 - fi = 3^11 + gj = 3 
for {pp\{db)v,A\B^) , where j = 2, ■ ■ ■ , 5 and g\\, gj^ and /|^ are expressed by 



(13) 



9\\ 



t3 



t3 



/J 



(14) 



We note that the power expansions of gi, fj ocl/t^ are due to the need of three hard gluons 



in the processes 



2J]. Two of them are used to create valence quark pairs just like the time- 



like baryonic form factors, and the third one is responsible for kicking the spectator quark 



in the B meson 



22(1 . We emphasize that the assumption of the only t-dependence for and 



fi in Eqs. (1121) . (fT3l) and ((Hj) may not be correct. In general, they are also functions of the 
other independent variable, such as {pp + puvY (Pp + PMvY- To verify our assumption, 
it is important to measure the baryonic leptonic B decays, such as B^ ppl^vi 22 1. 
From Eq. ([2]), the squared amplitude for B ppMy after summing over all spins is 
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given by 

= ^\aMy\'fM,\^[2ml{AE,Ep - t - m^^J + 2ml^X2t - SmJ + ml^^)]gl 
+24:171 M^mpniBiEp - Ep)gig2 + 2mp[{m% - tf - 2m\[^{m% + t)\gigi 
+AmBm.p{m\ - t - 2mljJ{Ep - Ep)gig^ 

+2m^^ [2m\{AEpEp - t - mlj^) - mj,^^(4mj - t - 2m\^^)\gl + 4ml^^mBt{Ep - Ep)g2g4 

+2m^^(t - 4m2)(m| - t - 2mlj^,)g2g^ + ^t[(m| - tf - Am\,^t\gl 

+2mB{ml - t)t{Ep - Ep)g4g^ + 2t[m|(Ep - Epf + mli^{t - 4ml)]gl 

+2[mi,^ + mliA< + + 2t) + m|(4i?pi?p- - t)]/^ 

+12ml^^mp{ml - t - 2m\j^)fif2 + AmBmp{m\ - t){Ep - Ep)fif4 

+8mp [mliEp - Epf + (t - Ami)] 

+2m^^ [2ml{AEpEp -t-mlj^)+ m\,^ (Sml + t + 2mU] fl 
+Am\j^mBt{Ep - Ep)f2U + 2m^,^^(t - 4mp)(m| - t - 2m\j^)f2h 
+\{t-Aml)[{ml-tf-Aml,^t\fl 

+2ms(t - Aml){ml - t){Ep - Ep)Uh + 2m^(t - AmlfEp - Epf fl^ , 



(15) 

where t = {pp + Ppf, Ep(^p) is the energy of the proton (antiproton). We note that in Eq. 
(IT^ . the terms related to C{B ppMy) have been ignored due to their small contributions 
to the branching ratios (< 1%). Replacing the proton and anti-proton energies -E'p(p) by the 
angle of 6 between the three-momenta of the vector meson and the proton in the dibaryon 
rest frame, i.e., 

mi, + t — m'i± cos 6 , ^ 

= s^T^ 

where /?p = (1 — Aml/tf/"^ and = -|- m\j^ + — 2m\j^t — 2m^t — 2m\j^m\, we can 
rewrite Eq. ( fT5l) as 

\A\^ = -f\aMy\^fM^{po + Pecose + pe2Cos^e), (17) 
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where 



Po 



m 



■Ah 



+ 16mjm^,^^ - 2(m| - 3t)m^^^, + (m| 



-12m 



+8mpmlj^{t - 4mJ)/i/5 



+ ^2 _ 2mlt]{gl + fl) - '^^^l^Mydl 

^ ^ ^ 2mi,^(m| + t) + (m| - tfMt + (t - 4m;)/2] 
+2mi,^(t - Aml)[tgl + {t - AmDfl] , 



2^ 



PpXy^ {2mij^^[6mpgig2 + t(fi'2fl'4 + /2/4)] 
-("^Mv -^B+ t)[2mp{gig5 + /1/4 - 2mpfif^) + ^(5(45(5 + /4/5)]} 



(18) 



The decay width F of -B — ppMy is given by 30 1 



dcos6 



dcos9 



+1 AmB-mMy)'^ 



1 J4m2 



COS^ . 



(19) 



(STrms)^ 

From Eq. f|T9|) . we can study the partial decay width dV /dcosO as a function of cos6', z.e., 
the angular distribution and define the angular asymmetry by 2J, |30|] 

T^d cose - r° -^dcos^ 

JO a cos & J — 1 CI cos 61 



Ae(Mv 



d cos ^ 

^1 dF 



rfT rfcOS^ + [° ;v^rfcOS^ ' 
JO a cos y J — 1 a cos 

which is equal to (A^+ — N^)/{N^ + A^_), where N± are the events with cos 6* > and 
cos^ < 0, respectively. From Eqs. ffT7|) . ffTUl) and (12U]) . we obtain 

I r{mB-iriMyf 

Ae{Mv) = - pedt. 



(20) 



(21) 



We can also define the direct CP asymmetry in i? ppMy by 25 1 

T{B- ppMy) - V{B+ ppM+] 



AcpiMv) 



T{B- ppMy) + T{B+ ppM^ 



(22) 



From Eqs. (fTSi) and (fT9i) . we get 



Acp{Mv) = 



\0lMv 


2 _ 




2 


\o(Mv 


2 + 


OiMv 


2 



(23) 



where denotes the value of the corresponding antiparticle. 



III. NUMERICAL ANALYSIS 



be 



n our numerical analysis, the vector-meson decay constants of fp and fx* are taken to 
311 



(/p,/i^*) = (0.205, 0.217) GeV, 



(24) 



and the form factors in the B My transition as functions of t are given by 32 1 

V^i[^](0) 



Vi[Ao]{t) 



(1 - t/M2p])(l - ait/Mj^j + a.tyMy ' 

^1,2(0) 



l-ait/M| + cT2tVM|' 
where the values of ^(0), Aq, 1,2(0), cri,2 and Mi, 2 are shown in Table [T] 0. 



(25) 

From the 



TABLE L Form factors of B ^ K* (p) at t = in Ref \32] with {Mi^Nh) = (5.37,5.42) and 
(5.27, 5.32) GeV for K* and p, respectively. 



B^K* (p) 


Vl 


Ao 


Ai 


A2 


f(0) 


0.44 (0.31) 


0.45 (0.30) 


0.36 (0.26) 


0.32 (0.24) 


0"! 


0.45 (0.59) 


0.46 (0.54) 


0.64 (0.73) 


1.23 (1.40) 


0-2 






0.36 (0.10) 


0.38 (0.50) 



experimental data in Refs. [l|, I33I, w e obtain the time-like form factor related coefficients of 
C\\ , Cj^ and D\\ in Eq. (ED to be [341 



C\i = 93AtliGeV'' , Cjr = -203.8+^^^;;; GeV^* , D\\ = 189. 91];^ GeV" . 



+103.6 



1+13.1 



(26) 



Here, we have assumed that C|| -jj and D\\ are real since their imaginary parts are expected to 



be small 



By fitting the measured branching ratios of i? ppK^*^ and 5° — ^ ppK*^ 
and angular distribution of i?^ ppK' 

-mm 

, the values of and in Eq. dH]) 
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are given by 24 1 

Nu = 124.0 ± 26.5 Gel/^ N, 



-200.2 ± 51.1 Gel/^ M,? 



-247.0 ± 48.7 Gel/^ 



M, 



-22.1 ± 15.7Ge■^/^ M,^ = 229.4 ± 22.7 Ge\/^ 



(27) 



Note that there are no terms related to M| ( or and /a) due to the relation of e -pMy = 0- 
To compare the form factors /j and Qi in Eqs. f|T2l) and f|T3l) . we show their central values 
in Table [TTl respectively. It is interesting to note that 



fi{B- 

MB' 
9.{B' 



pp) ~ -1.5/i(5° ->pp) , 

pp) ~ 1.55(1(5° ^ |)^) , 

pp) ~ 4/,(5°^pp), (2 = 2,4,5), 

pp) Ag,{B^-^pp), (i = 2,4,5) 



(28) 



In addition, we have 



/2 ~ -h > (<)/4, 



(29) 



for the charged (neutral) B transition. 



TABLE II: Central values of fi{gi) = t^Mgi) (in units ofGeV^). 





/i(5i)/"ip 


hih) 


74(54) 




B pp 
B^ pp 


149.3(291.4) 
-98.2(186.4) 


329.3(-329.3) 
-82.3,(82.3) 


29.5(-29.5) 
-7.3(7.3) 


-305.9(305.9) 
76.5(-76.5) 



In Figs. [T^ and[D3, we show the differential branching ratios of dBr/drripp and dBr / dcosO 
as functions of mpp and cos9 which demonstrate the apparent threshold enhancement and 
angular asymmetry, respectively. The inetrgrated results for the decay branching ratios and 
the angular and direct CP asymmetries are presented in Table IIIII 

Here we have averaged the particle and antiparticle contributions for the branching ratios 
and angular asymmetries. We note that in the table we have only included the errors from 
the experimental data of the three-body baryonic B decays, while the others such as those 
from nonfactorizable effects and CKM matrix elements can be referred to Ref. 1251. In 
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2 2.5 3 3.5 4 -1 -0.5 0.5 1 

mpp (GeV) oos9 



FIG. 1: Differential branching fractions of (a) dBr / dnipp and (b) dBr / dcosO as functions of rUpp 
and cosO, where the sohd, dot and dash hnes represent the decays of ppK*~ , B^ — > ppK*^ 

and B~ — > ppp~ , respectively. 

TABLE III: Decay branching ratios (in units of 10~^ j and angular and direct CP asymmetries 
in B — > ppMy (My = K*,p), where the errors are from the experimental data of the three-body 
baryonic B decays. 



My 


Br 


Ae{Mv) 


Acp{Mv) 




6.0 ± 1.3 


0.13 ±0.05 


0.22 ±0.01 




0.9 ±0.3 


-0.27 ±0.06 


0.013 ±0.001 




28.8 ±2.1 


0.11 ±0.06 


-0.029 ± 0.009 



addition, we have ignored the imaginary parts of the form factors in our determinations of 
the direct CP asymmetries. 

From Table IIIIl we see that Br{B^ ppK*^) = 0.9 ± 0.3 is consistent with the experi- 
mental value in Eq. ([T]). In addition, we can explain the inequality of Br{B~ ppK*~) > 
Br{B^ ppK*^) in terms of the QCD counting rules illustrated in Eq. fl28|) . which is 
also demonstrated in Fig. [1^. The large branching ratio of Br{B~ — ppp~) ~ 30 x 10~^ 
can be understood from the ratio of \VubV*^aifp/VthV^gaA,fK*\'^ ~ 4 as the baryonic B ^ pp 
transition form factors are the same between B~ ppp~ and B~ ppK*~ . It is in- 
teresting to point out that if Br{B~ ppK*~) ~ 10 x 10"^ as indicated by the central 
value of the Belle result in Eq. ([T]), we find that Br{B^ PPP ) ~ 50 x 10~^ while 
Br{B^ ppK*^) ~ 1.2 X 10"^ which is close to the central value of the BaBar data in 
Eq. ([1]). Since the main contributions to B^ ppp^ are color suppressed a2 — Oi/lOO 
with A'c = 3 as seen from Eqs. ([2]) and ([3]), the prediction of Br{B^ ppp^) is in the 
order of 10~^. Even if a2 ends up around the order of 10~^ with Nc = (2, oo) for possible 
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nonfactorizable effects, it merely shifts to be O(10~^). Although currently there are no data 
about the branching ratios involving p as well as the angular distributions, our predictions, 
in particular the large decay branching ratio for the charged p modes, provide good tests for 
the PQCD approach in the ongoing and future B factories. 

We remark that our results of A0{K*~p) ~ 10% and |ylg(i^*°)| ~ 30% are smaller than 
A0{K~) = 0.59;*^q;o7 measured by Belle j^. It is interesting to note that the large value of 
Acp{K*^) ~ 22% is in agreement with the BABAR data of (26 ± 19)% as given in Ref. fl- 
On the other hand, the direct CP asymmetries in ppK*^ and B —* ppp are too small 

to be measured. 



IV. CONCLUSIONS 



We have studied the three-body baryonic B decays of -B ^ ppMy in the standard 
model. The baryonic matrix elements have been calculated in terms of the SU{'i) flavor 
symmetry and the QCD power counting rules within the PQCD. Explicitly, we have 
found that the decay Brs, Aq and Acp of {B~ — > ppK*"^B^ ppK*^^B~ PPP ) 
are around (6,1,30) x 10"^ (13,-27,11)% and (22,1,-3)%, respectively. Our result of 
Br{B^ — > ppK*^) is consistent with the BaBar and Belle data, while that of Br{B~ — > ppp~) 
should be observed by the ongoing B experiments at the current B factories soon. Finally, 
we remark that our large prediction of Acp{K*'^) can be used to test the standard model 
and search for new physics. More precise measurements are clearly needed at the current 
and future B factories. 
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